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©f Importance in materials lifce net&ls at high temper®- 
tores* rubbers* plastics and most high polymers# 


Viscoelasticity theory allows for dissipation* I 
an offset which is not considered by classical elasticity 



solid to re wart to Its natural state* On the other hand 
a fluid cannot he in eeuilibyiw® under shear# It will 
flow if Stearin** stresses* hows war snail* are applied# j 

It does not hare a natural state and as such work done i 

on it by shear Is dissipated as heat# i 


Purely wi scows and purely elastic substance# ! 
are idealisations? theories based on these idealisations : 




*1 2 j 
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thermodynamic studies o<* til® foundation of the 11 .near 
theory; the re-examlna tl on of the linear theory ^mm 
the view point ©r non-'l incur continuum mech anlcs* The 
second category of investigations deals with methods 
of integration appropriate to the fundamental boundary 
value problems. An example being the correspondence 
principle which linlrs the linear theories of visco- 
elasticity and elasticity* So m@ physical chemists 
such as Ferry and fobolsky have been concerned with 
the linear response of viscoelastic materials with 
the aim of learning mm thing about the molecular 
constitution of matter* Bland and others have solved 
a large number of boundary value problems using diffe- 
rent formulations* mm general formulations have b mn 
proposed by mathematielans and those working in the 
fields of theoretical mechanics. 

•■#© will develop a linear theory of visco- 
elasticity which takes into account the effects of 

couple stresses* - . 
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i. s immk m s zM Eism 

l.l. „Qm, ,ai’S8asl,sna*,.-tecs3ls.t3iin. 

A generalised viscoelastic solid is spec If led 
fey the existence of a functional donation of state 
connecting stress < cr*5, strain (€}, time <t) and 
temperature CO) 

z^* - C < 3~? ^- ■> ^ > — O ,** 1,1*1 

Initially, for simplicity, the discussion is restricted 
to a solid stressed under uniaxial tension or compre- 
ssion* The appropriate generalisation to include three 
dimensional cases is given later* fh® relation 1.1*1 
■which may include time differentials and Integral opera- 
tors of arbitrary orders. Immediately excludes problems 
associated with the slip and dislocation of metals and 
fracture of solids* 

By expanding the functional relation 1*1*1 to 
the lowest order of non-vanishing terns, (linear in 
stress and strain) and by camming the temperature to 
enter only as a constant parameter! a relation for 
Isothermal behaviour can be derived Uhich conforms to 
the Boltsmann superposition principle for linear visco- 
elastic solids* 

The phyiieel phenorvna»i which distinguish 
viscoelastic substances from pure viscous fluids and 
pure elastic solids are recovery of strain or recoil 
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when the stress is released an# the flow or creep 
tinder cons taut stress respectively* 

Of the fssny pairs of associated stress- strain 
cycles* It is usual to choose either creep or relaxa- 
tion behaviour ss a fundamental measure the mechanical 
properties or & linear solid. Assuming the response to 
contain instantaneous components of stress and strain 
related hy s dynamic modulus the creep and relaxation 
h@hav5.onr may he expressed hy the linear relation 

£(4) = V'Ct-Tj^&ctr ... 1 . 1.8 

<rft? = £ ec-o ... 1.1.3 

-'V 

Sere y? it) and ^ < t) are the uniaxial creep and 
relaxation functions* 

and S? ( t) ate now non-aimensioimlised hy 
expressing 

y'C-^J - <^V/ 

*** 

*t*«re 8| is another instant having the dimension as 
that of 1 

In teres of non-dijsensionali sad creep lime- 
tion y>(t) and relaxation function ^ <t)« equations 
1*1*2 and. 1*1*3 are written down at * 

= ...l.i.# 

— % 

^ * * s Ci> - £,{*9 W 


3,* I*# 


— t 6 I* 


In writ inf down 1*1*5 and 1 * 1 . 6 , it Isas been assumed 
that B* 1 and are finite cuanti< les. While this 
&9smptSb& excludes some of file highly Idealised 
models ©f viscoelasticity (e»g* t Voigt solid, IswtOKiii® 
fluid). If is almost certainly correct for real solids* 

We will discuss at a later stage how these equations 
can he used with certain restrictions for Ifewtonlen 
fluids and Voigt Solids* In the following discussions 
these categories should not he Included unless specifi- 
cally mentioned* 



Fro® physical observations we note that during 
creep the strain increases at a constant value of stress, 
while during relaxation the stress decreases at a cons- 
tant value of strain* By considering equations 1*1*6 and 
1*1*6 and keeping in mind the types of physical behaviours 
described above, we conclude that t Cl) 9^ (t) wust be 
a none tonic increasing function of tine Cfewionii® fluid 
and %tgt Solid being exceptions) * 

(ii) Y {%) must be a monotonia function of 
tine (Wewtonian fluid being an exception}* 

(Hi) As <r <t) can never change sign, during 
relaxation, <p Ct) is bounded by its upper Ifnii^C^ ) 
where J Y<P ( 00 ) , K being equal to # 

^ Ct) and y-' Ct) ate related* This em be 
shown by faking the Replace transferee of equation* 
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1.1*5 and 1*1*6* Hie equation# obtained ai*® i 

<? C sj = ^ <*■* r cs) j ™ ... !.«.! 

a~cs) - £ e ~ ' s ^ ^Jjecsj ... j.B.g 

Bcrof valence ©f the above relations yield# 

/>- ^ ^ CS) J[ £'V^^ fcsjj = i ... 1 . 8.3 


which in the case of IS » reduces to 



Snuatlon 1.2.4 is due to Gross (2) and is of some 
Importance in ©©-relating the enter and relaxation 
behaviours of a typical viscoelastic material* 

Making ns® of the Abelian theorem (for 8 — MS t 
through positive vain®# of 8, JJ- s-^C-sj = <££ 

S-Z-O £ — >- ^ 

and using actuation 1*2*3 we get t 

^ 'Y'C^fJ^d 1 , 2*8 

Which in case of I » Ij reduces to i 
•// . — <p £ / ~f~ ys c°e>yj ■= •** 1 * 2*6 

From 1*2*8 we get ^ < ES^ L-~ *## l*t#f 

r '■'■ ^ /-/< 50 (foe; 

Coroeeting the ultimate values of the creep and relaxa- 
tion functions# For a partially relaxing solid % cp 0*0 <C 
'(perspex at ambient temperature) while for a completely 
relaxing solid (Folyisobutylcne at ambient temperature) . 
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XaSU, 

It Is of practical Interest to formulate tho 
typo of *? function tqnribolftvixiff different material 
behaviours* For solids which exhibit instantaneous 

4* 

elastic behaviour (Maxwell solid) f <p> <t) * (1- eT 5 

is the sinplost form for <p obeying all the restriction® 
discussed earlier* 


Substituting in ©nuatlon 1.1*6 we get 

-Ct-T) ]0 ie&h 

crOO ~ C-0 - PF X ( C / ~~ e ' ■ J jr l#^*i 

v~ @0 

It w® tab® €. <t) * e o m then at p * t It gives 
the Instantaneous response a® 0~ » 1 €: e * 

For t * Co it gives 

C T* C?°y ~ ^ ^-o ~~ e => ••* 1*3*2 

Complete relaxation takes place When I » Whereas 
the stress a® suites a constant value at t * on - in case 
of 1 > Eg* This has been shows in Fig* 1* 


a 

€ 



£j= ^opcQ 


•6 


U 

cr 



Fig* i 

Complete or Partial relaxation of Maxwoll ooll# 



It should he noted that In choosing pit) • 
(%m we have loosed a unitsue relaxation pattern 


to all Materials* To account for quick ’•©taxation or 
slow relaxation feehavlour the general form propose* 1st 

— Off 

, . Cpct) - _ ... 1.3.3 

0( s& €L/r>^ cl. nrr) cxrfe ttlcxC zz.czmsAci'y? t cohlcJx /& boSt/i^e. , 

hdju~&m$L* 


Consider the equation 1*1*5 defining strain 


at tin© t 


CC'O - & <5~C-0-t~ J_ 


1 



<b 


dr 1 




The general torn of \js (t) is ofetalmfi from the rela- 
tion given (in section 1.2} in the taplae® transform 
plane as _ _ y 

[k-£ i s<PM][£~+5< l 'r‘*J=i ...1.8.3 

where <£> (t) is given fey^Ct) *[/-& * (t ] ... 1.3.3 

The ©lass of solids know! as Voigt Soli at 
do not exhibit any relaxation phenowsna* Their instanta- 
neous response to suddenly applied load is nil, but 
they show a well defined creep phenomena. Equation 1.1*6 
cannot be used for such materials, whereas we can use . 
equation 1*1*6 with 1 ^ infinity* This gives the strain 
at any tine as i 

^ C^f) ~ £j_ J~ If? cf gr * * • 1.4*1 

Y'(t) may fee chosen keeping in mini the fact that 
such a solid possesses elastic and viscous elements in 
parallel in the "macroscopic level, and .as such creep 
can never fee infinite, fence the Creep function is 

“Cl ..Mi J|| " 

iWlflQWl,* . ' . , , 



Cr*®p ©f © typical Voigt Solid 


For elastic solid © <p (t) « 0 y y(t) * 0» For 
tills case ecpiatlons 1.1.5 and 1*1.6 raduc® to 
e^r£)= ^'a-CV 

... 1 . 0*1 


u LV" Jj v c ... 

Tha st&tafel* for* of 9 ? Ct) * § <t> 

*&•*§ S' (t) is tit® Birae*s Delta function W 

$ 

j r Jt - 1 t r&Mms ©©nation 1*5.2 


i 
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The negative slim can h# Interpreted by taking 1* 

x dt 

as the internal resistance acting opposite to the 
applied traction. It should he noted that <r(t j In 
equation 1.5*4 3s the shear stress and E^ has the 
dimension of viscosity 


The constitutive equations for the three 
dimensional case will now he considered. It will he 
assumed that the material is isotropic In behaviour. 
As an isotropic elastic solid is specified by two 
independent elastic const«|ts f the behaviour of an - 
isotropic viscoelastic material will be specified by 
two independent functions of time* There are many 
ways of doing this* Wo will take either the creep or 
relaxation functions as the basic functions. 

The constitutive equations proposed are t 

7f(t)=7\ %y e rr (V ~ \ ^ Jz 





W/W 

Jr 



•## 1 * 6#1 






*?£(&'>>+ - 2 /# 

/ - r _Z 1 / . / 


r„(i) *■ 


■M 




d2B<%/ r 

(s\+M,) 7 ‘* r 


c/KCP 


r J~dT~ C,/r ... 1.«.8 

~<b 
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Where Sjj Is the Kroneeker delta* p 0 * 9 ? i af ® 
relaxation functions ana y^_ are creep functions* 

Boltzmann (3) suggested that there is no elastic 
after-effect for purely dllational strains. This lilies 
that the relaxation ana creep functions are zero. Folsky 
ana sail. (4) proposed that this is probably a realistic 
assumption for a crystalline solid for which hydro- 
static pressures merely change interatomic distances. 
Another assumption discussed by Folsky and Shi la 
that and % * where 9? $ yj£ art relaxation 

and creep functions under hydrostatic condition end 9? t 
ijg mm those under pure shear respectively. This 
hypothesis is clearly false for Newtonian fluids under 
shear. 
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Si. 

Consider a body g> 0 bounded by the surface 38 e 
(as Sheen in Fig* 2*1)* Let the body <8 0 b® sobleeted 
to a system of forces, couple* s and displacements * I Sow 
Imagine a portion Bj of B 0 be separated from the rest of 
3 0 by the surface ''2*4 * Then, according to the generalised 

fox® of the stress principle 
of Euler and Conchy, the 
effect of portion on 

Sj may be represented fey a 
str®ss«>v@©t©r C%J and 
a eouple*st re ss«veeto r 'WJC'*) 
acting over the surface '^>Q i 
bounding S £ * Call any pair of distributions (jf, w ) and 
C , n? y ) of force-stress and couple-stress equipollent 

If they yield the son* mines for the resultant for©# 
mi moment ©n ©eery part of the material body S a • then* 

If the distribution ( t*”?) ©n surface ^4 be equipollent 
to the ©Mtaet action (due to distributions ( j£,w/ ) } of 

( 4>~4 ) « 4 » ^ li equivalent to saying that the 
part & 4 eaanet *fl»el* the difference then the portion 
( B d -B i > ii removed* ttstil now, most of the ee*tc in 
eontlnm mechanics has dealt with the son-polar ease* 
flie noswfolar case is characterised by the conditions 
W t 23 3 * 0 and ^ * © where is the body moment 
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«# 


per unit mass analogous to the body force per unit mass b * 
the basic principle ©f balance of moment of taomenitat with 


such a characterisation leads to the conclusion that the 

obtained from. 


stress tensor /,-• 

<J 


4-=^ 5- 1* .jra-otrlc 


1*®* ^ - 7]j * %en maple-stresses and body nonents 
are retained # the balance principle for angular momentum 


does not lead in general to the symmetry of stress tensor* 


In fact there is no reason for neglecting the 
couple stress rector on * There may be cases where Mjj 
obtained from on. - M.- ^ may haw® an insignificant 
magnitude or In certain classes of material® couple stress 
may not occur* 

Wmm Cauchy* s first law ? (balance of linear 
mmmtm) <«. <*a .« that tea «toM« t»n«r ^ «Ut> 
for all smte rials* Bat from Cauchy* s second law it cm 
be seen that ewen if the body eouple l exists# the field 
equations can be satisfied by taking couple stresses to 
be ssro* Urns Hie concept of body moment is essentially 
unrelated to the definition of couple* stress# fast ss the 
concept of body force is unrelated to the definition of 
force-stress. 

According to Truesdell and fnapln <6># body 
momenta mm first considered by Wawell tWTB}# barasor 
<1892)# Combebia© (1902) and Batten <H04)« 



•I 15 «» 


Couple stresses have Men considered by VolgtC 1887 ) , 
Em and F# Cosserai (1909), Heim (1913) and Gunther (1958)* 

Trnesdell and Toupin haw given & isodem fieri va* 
tlon of the Co as® rat equations, foupin (5) ■was the first 
to give a constitutive equation for couple stresses in 
elastic materials subject to finite deformations* Mlndlin 
and Tiers ten (7) linearised Tenpin * s constitutive equations 
and applied it to a number of problems such as effects of 
couple stresses cm vibration, wave motion and stress concen- 
tration* They haw also discussed the uniqueness of 
solutions* In a recent paper Stole®® (8) has proposed a 
theory for couple stresses 1m fluids* 

’#§ have made an attempt her® to extend the 
theory to viscoelastic materials* 

Eti i~:fi MmmIMJLjmaMmm 

The equations for the 
conservation of mass, the 
Balance of aM 

moment of tsonentna and 
the Conservation of 

*•** s#i*i 

jL-.[ V f dv - f k * J ~ ^ **• 2*1*8 

dt-J s r r J Sj J w 



A f fdv = a 
dt Ja 



#*$ 2JS I m 



where is the derivative , with aspect to tl mm* 
following a material particle* p is th© stass density* 
r is th© spatial position vector from a ftacsd origin, 
nj Is the materiai velocity , £ is th® ©sow 
per unit mass, 4 Is th® influx of energy per unit area 
and is th© energy source density per unit mass* 

In ©aeh of th© equations 2*1*1 to 2*1*4 the 
left hand tern gives th® quantity In balance, the 
surface integral In ih® right giro® the influx relative 
to the materiel and the volume integral gives It* s mam* 

Consideration of the equilibrium of forces 

,>w.. ■dktu.Jkm eaMSh* ^auJaieL » irifnnYrrfa ,m fi'irliftti- ■*»« fA|- ini m afihrfnwi rnni mm jwfcfa Jik. ma, .AseeaBa^., 1 ”^ eb, j ie^ w a&..^iia t .. sM k 

fitiQfijJIB OR fill ®A^8#TOfi¥^ | fill So# dvXVMi ojr 

the tetrahedron shrinks to aero, leads to the definition 
of the stress tensor ^ as 

~ **• 2*1*S 

in analogous treatment of the equilibria* of ament* 
noting on the tetrahedron yield* the definition of the 





**• 2 * 1*8 
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couple stress tensor M } . 

eJ 


m 


Equation 2*1,1 on simplification (see Appendix t) yield* 
the continuity squat ten 


* + 

f+ P u l ,1 ...2.1.7 

Hep® 2 A is tit® material velocity -vector# 

£ ‘ 

Using 2,1*5 and 2.1,7, equation 2*1.2 can be 
simplified to give the equation for balance of linear 
momentum ( Cauchy* s first lav of motion, see Appendix t) 
as s 

7** * ~t~ ~ #•* 2 * 1 *® 

J‘>J 

where Cl. is the material acceleration vector# 

♦file equation for the balance of angular mmmu* 
turn (Cauchy* s second lav of motion^ can now he written 
down, after simplifying equation 2*1*3 with the help of 

2,1.6 to 2*1*8 (see Appendix l) t as t 

€y K 7jK + M j:,j + = ° ... 3.1.B 

it can he seen from 2*1*8 that given an acce- 
leration field and a body force distribution, the stress 
tensor has to exist* On toe other hand equation 2*1*8 
can he • satisfied by talcing couple stress tenser © 
even if the body moment l* Is not sere* 
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Ths energy ©ma&ticm giw in conation 2*1*4 
on amplification toy using 2,1,5 to 2,1,9 (Appendix t) 
nlms t 


p£ ■= 7/f + M °j K u - ^ 


... 2 , 1,10 


yfh.®m ~Jj< is th© synjet-pie part of the stress tensor. 


,<Z> 


M;j U th. Oeviatortc part of oouplo *«« 


tensor given by 


M,'i ~ {f^'/ J Myr^Jj) *** 


1 1 % 


sZ> 

y = ^ 

It is interesting to not# that only th# symmetric part 
©f stress tensor and the deiriatorie part of soupls stress 
tensor appear in the energy equation. 


be. 


f sating th® curl of equation 2,1,9 m git 

V* [e ;Jk Tj K+ M jy +?(;] 

(’—mm ^ fj K J 

fa„- r ^),„ = e ^' %■» * 


£L, e iA 17*>” -f <?«-»/ M Ji’j» -+ e mr,:^)’» 0 


-T- A 
'vrY7,7J 


Sines 




**• 


7X F. M - VA V • M 


D 


Cauchy* 3 first 1m of motion 2 . 1*8 12000063 * 


fa, = Tptj +£[' -.jk 7 - * e lJK (P4 4 - j + f4 


y 


### 2#X#!U3- 


vh«™ Tj* = iCr J; + Zy) 


% !e can also pit the antisymmetric part of the stress 
tensor from conation 2 * 1 . 9 * Urns multiplying through 


fey £L„; «® git e 




e y'/«c Z'K 


-h G 


yyml 


/V-‘ < 

Jh J 




T^ n = - -kfemnl M jy + ^S! *** Sa * 14 

If lit fee the components of displacement ©? a 
material particle f then for infinitesimal displacements 9 

Xfia O011f®I»OTS Or SCmiB &TO glTBfi DJr I 

€;• = / (f^y ^ ... 8*2.1 

fli# components of the strain tensor glron fey this smell 
st rein agnation for * ? */ are half the physical shear 
strain co®pon©Bt» • 

the material rotation rector is defined 


**s SO i 


The three components of 03- give the snail 
rotation of the material particle about the three 
ares. We define the cursnaiure-twl st tensor K,y as 
the gradient of rotation vector Ca)j 

i#e* ~ ~ . 2 " si 00 * 2*8.3 

Three diagonal components of K !f give the 
twist of the material per unit length about the three 
exes* On the other hand, the off-diagonal components 
of A C)j are a measure of the curvatures indued in 

the various planes# 


The rates of strain t*nsor y the vertieiiy and 
the curvature twist rate tenser are given respectively 
by * 


~ i ( Ul >j * 

»*# 2*2*4 

co/ - 2 e i/ x ^k,j 

• ** 2*2*0 

K-;j = - i e y^ ^5/ 

I*** 2*8*0 

from 8*2*3 ani 2*2.0 that 

AT// * o , >£// = o 

**# 8#8#T 


Defining tiie antisymmetric part of the veloeiiy 


gradient £ C ~ the vertieity teaser* 

it mm he en tfoefc t 


W/j ~ ! ^y /c 

*k = i ^ % 


*** 


2 . 2.8 


@##*0 




MSn&ln and Tiersten ( T ) h&m g iwn 
linear constitutive equations for the polar cast in 
elasticity* Stoke* C Q ) has given tho gnu© for 
fluids* W# will um these results to propose mm«* 
titutiee equations for the slseoelettle ease* 

Considering the esastitutlee equations prop osed 
ky Kiadlin sad fiarstoa and Stokes ws propose the 
following ©oastituti'TO equations for the synraeirte 
part of stress teaser and dorlatoric part of couple 
stress tensor for isotropic elseoslastio materials i» 

77/ =-Pie ry %. -t-iM €;■ - e r cz) 4 T 

. 71171 . , . , 

-ZMjJ *?, &-T) 6,y (T) cl? 


II *#\ s#i*i 
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functions 9^ , CP, cp z> cp 3 are re luxation 
functions with fame similar to those giver* in Section !• 


The constants A 5 A ? $ % y % 


/ t of t P7 7 f may be constants or in general 


fZdtW of 49 etc. 

P\ ? ?/fj have dimensions of Force/ Area !•« 


/ have dimensions of Force l«e* 





I’l ’%'% 1 ?* 

The ratios between the second set and first set of 

§ 

constant® have the dissensions of V * 


The antisyssaetrie part of stress tensor and 
the trace of couple stress tensor are loft tmdetesoiised 
by the constitutive equations and can be detonstosi 
frws the b©t»da?y conditions# m tmu® already seen 
that thane Quantities do not ewe in In the energy 


of all the eduations 


far we haw coat across the fallowing 


equations t 

f + = O 

Ey * ?h = PO; 


one equation 


n*u? 


Three equations 2,1 # S 


**f 23 $ 


^ 'J ^ ~0k Mji 9 j Thm® equations 


*## 2*1*0 


£ - 
7 


S' 


S £y- i -tf J .Ky-h / ,;+f$ 

One equation 

Z ( tL: j + U j>0 

Six Independent ©qns 

i C*:,j+«j.iJ 

Six independent eqns, 


Ki;~ 


' J J 9C Z d sir t»Si Hint aquations 


• #* 2.1.1G 

» 2*2*1 

► 2*2*0 

2*2*3 


* » ~ —-C. Zl 

‘J * £ j st- t-,si fin© equations *** 2 * 2*0 
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